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Spin-dependent transport through a multilevel quantum dot weakly coupled to ferromagnetic leads 
is analyzed theoretically by means of the real-time diagrammatic technique. Both the sequential and 
cotunneling processes are taken into account, which makes the results on tunnel magnetoresistance 
(TMR) and shot noise applicable in the whole range of relevant bias and gate voltages. Suppression 
of the TMR due to inelastic cotunneling and super-Poissonian shot noise have been found in some of 
the Coulomb blockade regions. Furthermore, in the Coulomb blockade regime there is an additional 
contribution to the noise due to bunching of cotunneling processes involving the spin-majority 
electrons. On the other hand, in the sequential tunneling regime TMR oscillates with the bias 
voltage, while the current noise is generally sub-Poissonian. 
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Introduction: Transport properties of quantum 
dots coupled to ferromagnetic leads are currently 
a subject of extensive experimental and theoretical 
studiesi^^^^^i^^aaiiO^asOi^ This interest is stim- 
ulated by expected applications in spintronics and quan- 
tum computing. ^'^'^ When a quantum dot is coupled to 
two ferromagnetic leads, its transport properties depend 
on the magnetic configuration of the system. This is the 
so-called tunnel magnetoresistance (TMR) effect, which 
is characterized by the ratio TMR = (/p — /ap)/^ap, 
where Ip (Zap) is the current fiowing through the system 
in the parallel (antiparallel) configuration. ^^^^^^ When 
the coupling between the dot and leads is strong, the 
Kondo physics emerges for T < Tk, where Tk is the 
Kondo temperaturei^ In turn, single-electron charging in 
the weak coupling regime leads to the Coulomb blockade 
phenomena. Sequential (first order) transport is expo- 
nentially suppressed in the blockade regime. The current 
fiows then due to cotunneling (second order) processes 
involving correlated tunneling through virtual states of 
the dot, whereas outside the blockade regime transport 
is dominated by sequential tunneling processes. So far 
fully systematic considerations (taking into account both 
sequential and cotunneling processes) of spin-dependent 
transport through a quantum dot in the weak coupling 
regime have been restricted mainly to single-level quan- 
tum dots.^ In real systems, however, usually more than 
one energy level participate in transport, leading to more 
complex and interesting transport characteristics J^^i^ It 
has been shown recently that the Fano factor in the 
Coulomb blockade regime calculated in the first order ap- 
proximation is larger than unity Moreover, the TMR 
was found then to be independent of the gate voltage. In 
this paper we extend the theoretical studies by including 
the cotunneling processes, and show that the shot noise 
in the Coulomb blockade regime can be super-Poissonian, 
although the Fano factor is significantly reduced by the 
cotunneling processes. Apart from this, we show that 
the TMR in the blockade regime is considerably modi- 



fied by cotunneling processes, and can be either enhanced 
or reduced in comparison to that in the first order ap- 
proximation, depending on the transport regime. Our 
considerations are based on the real-time diagrammatic 
technique^'^^^^^ which, after taking into account both the 
first and second-order contributions, allows us to ana- 
lyze transport in the full weak coupling regime, i.e. in 
the cotunneling, cotunneling-assisted sequential, and se- 
quential tunneling regimes. Furthermore, to calculate 
the shot noise in the cotunneling regime, we include the 
non-Markovian effects, which were neglected in previ- 
ous considerations.^^ 

Model: We consider a two-level quantum dot weakly 
coupled to external ferromagnetic leads whose magnetic 
moments are either parallel or antiparallel. The Hamilto- 
nian of the system reads, H = Hi, + i^R + Hp, + Ht- The 
first two terms describe noninteracting itinerant electrons 
in the leads, Hr = Xlkcr ^n^acl^^^Cri^a for the left (r = L) 
and right (r = R) leads, where cl^^ (crkcr) creates (an- 
nihilates) an electron with the wave vector k and spin a 
in the lead r, and e^kcr is the corresponding dispersion 
relation. The quantum dot is described by 



^ ^ ^jT^'i + X] ^la^2a' , (1) 



where rija = dj^^ja and dj^ {dja) is the creation (anni- 
hilation) operator of an electron with spin cr in the jth 
level (j = 1,2), £j is the corresponding single-particle 
energy, and U (U') is the on- level (inter- level) Coulomb 
repulsion parameter. The tunnel Hamiltonian, takes 



the form: H^ = Er=L,RE\.jaitrjCl^Jja 



where trj is the relevant tunneling matrix element. Cou- 
pling of the jth level to the spin- majority (spin- minority) 
electron band of the lead r is described by = 



27r|t 



r-)A 



I'^pt'^ ^ = r^_^-(l ±Pr), where F^^^ , ^^^^ 

while pt'' ^ and Pr are the spin-dependent density of 
states and spin polarization in the lead r, respectively. 
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In the following we assume Trj =T/2 and Ph = Pr = P- 
Method: In order to calculate the spin-polarized trans- 
port through a two-level quantum dot in the sequential 
and cotunneling regimes, we employ the real-time dia- 
grammatic technique, ^'^^'^^ which consists in a system- 
atic expansion of the quantum dot (reduced) density ma- 
trix and the current operator with respect to the dot-lead 
coupling strength F. The current operator / is defined as 
/ = (/r-/l)/2, with ir = -i{e/h) T.j{^rjcl^Jja- 

t^jdj^Cr\^a) being the current fiowing from the dot to the 
lead r. Time evolution of the reduced density matrix 
can be visualized as a sequence of irreducible self-energy 
blocks, on the Keldysh contour. The matrix ele- 

ments W^^f describe transitions between the many-body 
states Ix) and \x') of the two- level dot.^^ The full prop- 
agation of the dot density matrix is given by the Dyson 
equation, which is further transformed into a general ki- 
netic equation for the elements of the reduced density 
matrix. With the aid of the matrix notation introduced 
in Ref. [21], all the quantities of interest can be defined 
in terms of the following self-energy matrices: W, W\ 
W^^, 9W, and dW^. The matrix W^*^^^^ is the self-energy 
matrix with one internal vertex (two internal vertices) re- 
sulting from the expansion of the tunneling Hamiltonian 
replaced by the current operator, while dW and dW^ 
are partial derivatives of W and with respect to the 
convergence factor of the Laplace transform. Using the 
above matrices, the stationary occupation probabilities 
can be found from, (Wp^*)^ = TS^^^ , where p^* is the 
vector containing probabilities and the matrix W is given 
by W with one arbitrary row xo replaced by (F, . . . , F) 
due to the normalization, J^x^x ~ "^^^ current fiow- 
ing through the system can be then found from 

/=^Tr{WV*}- (2) 

Finally, the zero- frequency current noise, S = 
2 dtamm + mm - 2{i)^), is given by^i 

S = —Tv { [W" + (PW' + p'* <E> e'^dW^)] p'*} , 

(3) 

where the object P is calculated from: WP = l(p^* (g) 
— 1 — 9Wp®* e^), with 1 being the unit vector with 
row Xo set to zero, and = (1, . . . , l).^"*^ 

To calculate the transport properties order by order in 
tunneling processes, we expand the self-energy matrices, 
^(i,ii) ^ w(i'")^i)+W(i'")^2)^ ^^^^ occupations, 

pst ^ pSt(O) ^ pst(l) ^ ^ ^ ^ ^ ^^^^ p ^ p(-l) ^ p(0) ^ ^ ^ ^ ^ 

respectively. The self-energies can be calculated using 
the corresponding diagrammatic rulesj^*^ The first or- 
der of expansion corresponds to the sequential tunneling, 
whereas the second one to cotunneling. Thus, by taking 
into account all the first and second-order contributions, 
we are able to resolve transport properties in the full 
range of the bias and gate voltages. 

Results on TMR: The TMR as a function of the 
bias voltage and position of the dot levels is shown 



in Fig. [H^b). To facilitate the identification of differ- 
ent transport regimes, we present in Fig. [Tfa) the den- 
sity plot of the differential conductance G in the par- 
allel configuration (differential conductance in the an- 
tiparallel configuration is qualitatively similar^^). Since 
position of the dot levels can be shifted with a gate 
voltage. Figs. [Tfa) and (b) can be viewed as a bias 
and gate voltage dependence of TMR and G, respec- 
tively. By sweeping the gate voltage in the linear re- 
sponse regime, the number of electrons in the dot can 
be changed successively. More precisely, this happens 
when: s = ei = 0, s = -{5e + [/'), e = -{U + U'), and 
e = —{Se-\-U + 2/7'), where Se = €2 — Si is the level 
spacing. Thus, for £ > [-{Se + /7 + 2/7') > e], the dot 
is empty (fully occupied). When > £ > — (fe + U') 
[-{U ^ U') > e > -{Se ^ U ^ 2U% there is a single 
electron (three electrons) in the dot. On the other hand, 
for -(fe + /7') >£>-([/ + [/'), the dot is occupied by 
two electrons, one on each orbital level. 

In the case of empty (fully occupied) dot in the 
Coulomb blockade regime the current fiows only due to 
elastic non-spin-flip cotunneling processes. Such pro- 
cesses are fully coherent and do not affect the charge and 
spin state of the dot. As a result, the system behaves as 
a single ferromagnetic tunnel junction, yielding the TMR 
given by the Julliere formula,^^^^ TMR = 2pV(l - p'^). 
However, in the other blockade regions both the non- 
spin-flip and spin-flip cotunneling processes are allowed, 
leading to the suppression of the TMR. This can be seen 
in Fig. [TJc), where we plot the gate voltage dependence 
of the TMR in the linear response regime. When there 
is a single electron on one of the two orbital levels, the 
TMR is decreased. This is due to the spin-flip cotun- 
neling processes which provide a channel for spin relax- 
ation, decreasing the difference between conductance in 
the parallel and antiparallel conflgurations and thus re- 
ducing the TMR.^ On the other hand, in the regime when 
each dot level is singly occupied (doubly occupied dot), 
the amount of spin-flip cotunneling is increased and the 
suppression of TMR is even more pronounced. However, 
when the bias voltage is increased, the central minimum 
in TMR [Fig.[TJc)] transforms into a local maximum due 
to a nonequilibrium spin accumulation in the dot, see 
Fig.ffl^d) and (e). 

The bias voltage dependence of the TMR for several 
values of the level position is displayed in Fig. [T](f)-(h). 
When the dot is empty in equilibrium {e = 2 meV), the 
TMR in the cotunneling regime is given by the Julliere 
value. However, once the bias voltage reaches the thresh- 
old voltage {V ^ 3 mV), the sequential processes are 
allowed and TMR drops, see Fig. Hff). If the ground 
state is singly occupied {e = —1.75 meV), a nonequilib- 
rium spin accumulation in doublet states (p^^q) ^UO)) 
is built up with increasing bias voltage*^ This leads to 
an enhanced TMR, which again starts to drop around 
the threshold for sequential tunneling, see Fig. [T](g). Fi- 
nally, in Fig. [IJh) we show the bias dependence of TMR 
for the case when the dot is doubly occupied in equi- 
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FIG. 1: (Color online) The differential conductance G in the parallel configuration (a) and the TMR (b) as a function of the 
bias voltage V and level position e = £i. Parts (c)-(e) show the TMR as a function of e for several values of V, whereas parts 
(f)-(h) display the bias voltage dependence of TMR for several values of s. The dashed lines in (c)-(h) show the first-order 
contributions. The parameters are: /cbT = 0.15 meV, S2 — = 1.5 meV, [/ = 5 meV, U' = 2 meV, F = 0.1 meV, and p = 0.5. 



librium = —5.25 meV). Now, the spin accumulation 
in triplet states {pf^^^^ 7^ ^Ui)) ^i^e to an increase 

in TMR.^^ We also note that in the transport regime 
where the sequential tunneling is allowed, more and more 
charge states become active in transport with increasing 
the bias voltage. This gives rise to step-like I — V char- 
acteristics and the oscillatory-like behavior of the TMR, 
see Fig.[I{f)-(h). 

In Fig.[T](c)-(h) we also showed the TMR calculated in 
the first order (sequential) approximation (dashed lines). 
It is evident that the role of cotunneling processes is par- 
ticularly pronounced in the blockade regions, where the 
cotunneling processes dominate over the sequential ones, 
and lead to a significant enhancement (or reduction) of 
TMR. Outside the blockade regions, TMR is determined 
mainly by sequential transport, so the contribution from 
cotunneling processes is rather minor. 

Results on shot noise: Upon calculating the current / 
and the zero- frequency current noise 5*, one can deter- 
mine the Fano factor F, F = 5'/(2e|/|). The Fano factor 
describes the deviation of 5* from the Poissonian shot 
noise given by Sp = 2e\I\. The bias and gate voltage de- 
pendence of the Fano factor in the parallel and antipar- 
allel magnetic configurations is shown in Fig. [2](a) and 
(b), respectively. When \eV\ < /cb^, S is dominated by 
thermal Nyquist- Johnson noise, otherwise the noise due 
to the discrete nature of charge (shot noise) dominates 
In the limit of F ^ 0, the current tends to zero whereas 
the current noise is dominated by thermal noise. This 
leads to a divergency of the Fano factor in the linear re- 
sponse regime. Therefore, in the density plots of F we 
mark the low bias voltage regime with a black line. 

In the cotunneling regime, where the dot is empty (or 
fully occupied), the shot noise is Poissonian, Fig. [2](c). 



This is because the current flows then only due to elastic 
non-spin-flip second-order processes. ^^^^^ However, in the 
Coulomb blockade regions, where inelastic cotunneling 
processes also contribute, we find a pronounced super- 
Poissonian shot noise, see Fig. [SJ^d) and (e). This in- 
creased shot noise is related to bunching of electrons car- 
ried by different types of cotunneling processes. Further- 
more, in the case of magnetic leads there is an additional 
contribution to the noise coming from the difference be- 
tween the spin-up and spin-down channels which leads to 
bunching of fast cotunneling processes involving the ma- 
jority electrons. This is more pronounced in the parallel 
configuration where the difference between the two chan- 
nels is approximately equal to (1 +p)^/(l — p)^, while for 
the antiparallel configuration the two channels are com- 
parable, see Fig.[2](a) and (b). The variation of the Fano 
factor with the bias voltage is displayed in Fig. [21(f)- (h). 
When the dot is empty [Fig.[2](f)], the Fano factor is Pois- 
sonian and starts to drop at the onset of sequential tun- 
neling. With increasing bias voltage, the noise becomes 
super-Poissonian in some range of the transport voltage. 
In the case of singly and doubly occupied quantum dot, 
the bias voltage dependencies of F are qualitatively sim- 
ilar. In the Coulomb blockade regime the current noise 
becomes super-Poissonian (up to F 2.5 in the parallel 
configuration), indicating the existence and role of inelas- 
tic and spin-flip cotunneling processes. For example, in 
the case of doubly occupied dot. Fig. [21(h), the maximum 
Fano factor is found for voltages where inelastic cotunnel- 
ing between the two orbital levels is allowed, \eV\ ~ 2Se. 
On the other hand, in the sequential tunneling regime 
the Fano factor becomes suppressed and is generally sub- 
Poissonian due to the Coulomb correlations in sequential 
transport For comparison, in Fig. [2](c)-(h) we show 
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FIG. 2: (Color online) The Fano factor F as a function of the bias voltage and level position for the parallel (a) and antiparallel 
(b) magnetic configurations. Parts (c)-(h) show F in the parallel configuration as a function of the level position [(c)- (e)] and 
the bias voltage [(f)-(h)]. The dashed lines present the first-order contributions. The parameters are the same as in Fig. [1] 



the results calculated within the sequential tunneling ap- 
proximation (dashed lines), which clearly show the role 
of cotunneling processes in shot noise, particularly in 
the blockade regions. We also note that the occurrence 
of super-Poissonian shot noise in the Coulomb blockade 
regime has also been reported experimentally in quantum 
dots coupled to nonmagnetic leads. ^^^^^^^^ 

In conclusion, we have discussed the spin-polarized 
transport through a two-level quantum dot coupled to 
ferromagnetic leads in the sequential and cotunneling 
regimes. We have analyzed the dependence of the TMR 
and Fano factor on the bias and gate voltages. In the 
Coulomb blockade regime we have found a suppression 
of TMR for singly and doubly occupied quantum dot. 



Furthermore, in these transport regimes the inelastic co- 
tunneling processes lead to super-Poissonian shot noise, 
irrespective of magnetic configuration of the system. On 
the other hand, in the sequential tunneling regime, the 
current noise is generally sub-Poissonian, while TMR ex- 
hibits oscillatory-like dependence on the bias voltage. We 
also notice that generally the sequential tunneling ap- 
proximation in the Coulomb blockade regime underesti- 
mates the TMR and overestimates the Fano factor. 
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